Quantum Critical Behavior and Possible Triplet Superconductivity in Electron Doped 
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Density functional calculations are used to investigate the doping dependence of the electronic 
structure and magnetic properties in hexagonal Na^CoCh. The electronic structure is found to 
be highly two dimensional, even without accounting for the structural changes associated with 
hydration. At the local spin density approximation level, a weak itinerant ferromagnetic state is 
predicted for all doping levels in the range x = 0.3 to x = 0.7, with competing but weaker itinerant 
antiferromagnetic solutions. The Fermi surface, as expected, consists of simple rounded hexagonal 
cylinders, with additional small pockets depending on the c lattice parameter. Comparison with 
experiment implies substantial magnetic quantum fluctuations. Based on the Fermi surface size and 
the ferromagnetic tendency of this material, it is speculated that a triplet superconducting state 
analogous to that in S^RuCU may exist here. 



The last few years have seen the discovery of a num- 
ber of novel unconventional superconductors associated 
with magnetic phases. Focusing on triplet (or likely 
triplet) superconductors these include UGe2 (T c ~ IK), 1 
URhGe (T c ~ 0.25K), 2 and ZrZn 2 (T c ~ 0.3K), 3 ' 4 
where ferromagnetism coexists with superconductivity, 
and Sr2RuC>4 (T c ~ 1.5K), 6,7 which has a paramagnetic 
Fermi liquid normal state, but is "near" magnetic phases. 
Although the exact pairing mechanism has not been es- 
tablished in these materials, it is presumed that spin fluc- 
tuations are involved, most probably the quantum critical 
fluctuations in the materials with co-existing ferromag- 
netism and superconductivity. 8-14 In Sr2Ru04, strong 
nesting related antiferromagnetic spin fluctuations are 
found in local density approximation (LDA) calculations 
and experiment. 15 ' 16 In addition ferromagnetic fluctua- 
tions may also be present, and if so these would favor a 
triplet superconducting state. 17 

Recently, Takada and co-workers reported the discov- 
ery of superconductivity with T c ~ 5K, in the C0O2 
layer material, Na^CoO^ • yhbO. 18 Based on the two 
dimensional 3c? transition metal oxide structural motif, 
they speculate that the superconductivity may be re- 
lated to that of the cuprate high-T c superconductors. 
Here an alternate possibility is discussed, that is the 
connection with the triplet superconductors mentioned 
above. The parent compound of Na^CoO^ ■ 2/H2O is 
Nao.5Co02, which was synthesized and studied in single 
crystal form by Terasaki and co-workers. 19,20 The com- 
pound consists of C0O2 layers (Co is in at the center 
of distorted O octahedra formed by layers of O above 
and below the Co planes) separated by layers of of Na 
ions. 21 The result is a highly anisotropic metal with 
a large thermopower at room temperature and a sub- 
stantial renormalization of the Fermi liquid properties 
at low temperature. 19,20,22 Modifications of this com- 
pound consisting of substitutions for the Na layer have 
been made to obtain technologically useful thermoelec- 
tric materials, 23-26 and, in fact, Nao.sCoO^ is itself a 
good thermoelectric at high temperature. 27 Significantly, 



the Seebeck coefficient of Nao.sCoO^ as calculated from 
the LDA band structure using standard kinetic transport 
theory is in quantitative agreement with experiment, in- 
dicating that the Fermi liquid description at higher tem- 
perature is robust. 28 Band structure calculations for the 
higher performance so-called misfit compounds formed 
by replacement of the Na layer also reasonably account 
for room temperature experimental data. 29 

Turning to the low temperature properties, as 
mentioned, there are substantial renormalizations in 
Nao.sCoG^, for example in specific heat and suscepti- 
bility. LDA calculations indicate a ferromagnetic ground 
state for this parent compound, while experimentally fer- 
romagnetism is not observed, one possible explanation 
being that quantum fluctuations suppress the ordered fer- 
romagnetic state. 28 Interestingly, the related misfit com- 
pound, Ca3Co40g shows a large negative magnetoresis- 
tance at low temperature reaching 35% at 7T, which 
could be interpreted as a sign of proximity to a ferro- 
magnetic quantum critical point. Other evidence for a 
proximity to magnetism is the formation of a spin den- 
sity wave when Cu doped, 20 and the production of a weak 
magnetic ground state with strong magnetoresistance 
upon increasing the Na concentration to x = 0.75. 30 

The reported crystal structure of superconducting 
Na^CoO^ • 2/H2O differs from that of the parent com- 
pound by a reduction in the Na content to x ~ 0.35 and 
the intercalation of water molecules around the Na sites. 
The structure of the C0O2 sheets is very similar to that 
of the parent compound Nao.5Co02 - the in-plane a lat- 
tice parameter is slightly lower (2.823 A vs. 2.84 A) 
as is the height of the O (0.885 A vs. an LDA value 
of 0.908 A in the parent). It should be noted that a 
slight contraction of the C0O6 octahedra, as implied by 
these structural differences, is consistent with the lower 
Na content, in particular a lower electron count on the 
C0O2 sheets would be expected to lead to a small con- 
traction. 

Here, well converged LDA calculations are reported for 
Na x Co204 for £=0.3,0.5,0.7. In addition, calculations 
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are reported for a strained lattice corresponding to the 
structure reported for superconducting Na^OoC^ • J/H2O, 
but neglecting the intercalating water. The calculations 
were done using the general potential linearized aug- 
mented planewave method with local orbitals, 31 ' 32 as de- 
scribed in Ref. 28, except that better Brillouin zone sam- 
ples, corresponding to a minimum special k-points mesh 
of 16 x 16 x 2 in the hexagonal zone is used here to ob- 
tain convergence of the magnetic energies. As discussed 
in Ref. 28, a virtual crystal method is used to account 
for the partially occupied Na site. 

Although there is some hybridization, the valence band 
structure of Nao.sCo0 2 consists of three manifolds of 
bands separated by gaps - a lower lying occupied O 2p 
derived manifold, followed by Co t 2g and e g manifolds. 28 
As expected from ionic considerations the Fermi energy, 
Ep lies near the top of the t 2g manifold, which contains 
0.5 holes per Co ion. Because of the actual axial site 
symmetry, the t 2g manifold can be regarded as consist- 
ing of two-fold (also labeled e g ) and one- fold (labeled a g ) 
crystal field states. These overlap, but the top of the 
t 2g manifold is primarily of a g character in Na .5CoO 2 , 
with the result that the band structure near Ep can be 
roughly viewed as consisting of one band per Co ion with 
an filling of 3/4 (1/2 hole per Co). The present LDA cal- 
culations for the electronic structure of Nao.7Co02 and 
Nao.3Co02 follow this picture. These were done hold- 
ing the crystal structure fixed at that of Nao.sCoC^ and 
varying the Na occupancy in the virtual crystal. The 
dominant a g character at the top of the t 2g manifold is, 
however, lost for the strained lattice, discussed below. 

The magnetic properties are similar for the various Na 
contents in the range x = 0.3 to x = 0.7. Results of fixed 
spin moment constrained LDA calculations are shown in 
Fig. 1. In particular, itinerant ferromagnetism is found. 
In each case, the energy decreases with magnetization 
until a magnetization at which the band edge is reached 
in the majority spin. Then the energy increases rapidly 
reflecting the crystal field induced gap between the t 2g 
and e g manifolds. 

Thus, independent of x in this range, the LDA pre- 
dicts a ferromagnetic ground state, with a spin moment 
per Co equal to the number of holes (p = 1 — x) and a half 
metallic band structure (here we refer to the hole concen- 
tration as the concentration of holes in the t 2g manifold; 
without any Na p=l; Na electron dopes the sheets, which 
leads to a reduction in p) . The fixed spin moment curves 
show a shape crossing over from parabolic at low moment 
to more linear as the band edge is approached. As may 
be seen, the shapes and initial curvatures for the differ- 
ent doping levels are roughly similar. The trend towards 
slightly weaker initial curvatures at higher Na concentra- 
tion is possibly an artifact due to the fixed crystal struc- 
ture used in the present calculations. It reflects increas- 
ing hybridization (increasing band width and decreasing 
density of states) as charge is added to the C0O2 planes. 
In reality, the lattice would be expected to expand, per- 
haps compensating this trend. In any case, for this range 



of x, ferromagnetism with a magnetic energy of approx- 
imately, E(FM) f« 50p in meV/Co and spin moment 
M(FM) = p in pb/ Co is found. As mentioned, calcula- 
tions were also done for a strained cell with the structure 
of the superconducting sample, but without H 2 0. These 
calculations were done for £=0.5 and £=0.35, the latter 
corresponding to the experimentally determined doping 
level. In both cases, the LDA predicted a ferromagnetic 
state. The magnetic energies were E(FM) =20 and 27 
meV/Co for £=0.5 and £=0.35, respectively. Thus the 
behavior is similar, but the magnetic energies are some- 
what larger. 

LDA calculations were also done for an antifcrromag- 
netic configuration with the unit cell doubled along one 
of the in-plane lattice vectors. Thus, within a Co plane, 
each Co ion has four opposite spin nearest neighbors and 
two like spin nearest neighbors. At all three doping levels 
investigated an antiferromagnctic instability was found, 
but this instability is weaker than the ferromagnetic one. 
Details of the LDA moments and energies are given in 
Table I. Essentially, the energy of the antiferromagnetic 
configuration examined tracks the ferromagnetic energy 
at a value ~ 1/4 as large. 

The LDA generally provides a good description of itin- 
erant ferromagnetic materials. It is known to fail for 
strongly correlated oxides where on-site Coulomb (Hub- 
bard) repulsions play an important role in the physics. 
In such cases, the LDA underestimates the tendency of 
the material towards local moment formation and mag- 
netism. Here, the LDA is found to predict ferromagnetic 
ground states for materials that are paramagnetic met- 
als in experiment. While materials in which the LDA 
substantially overestimates the tendency towards mag- 
netism are rare, a number of such cases have been re- 
cently found. These are generically materials that are 
close to quantum critical points, and include Scaln, 33 
ZrZn 2 , 5 and Sr 3 Ru 2 7 (Ref. 34). Sr 3 Ru 2 7 displays a 
novel metamagnetic quantum critical point, 35 while, as 
mentioned, ZrZn 2 shows coexistence of ferromagnetism 
and superconductivity. 

Density functional theory is in principle an exact 
ground state theory. It should, therefore, correctly de- 
scribe the spin density of magnetic systems. However, 
common approximations to the exact density functional 
theory, such as the LDA, neglect Hubbard correlations 
beyond the mean field level, yielding the underestimated 
magnetic tendency of strongly Hubbard correlated sys- 
tems. Overestimates of magnetic tendencies, especially 
in the LDA are very much less common. Another type 
of correlations that is missed in these approximations are 
quantum spin fluctuations. This is because the LDA is 
parameterized based on electron gases with densities typ- 
ical for atoms and solids. However, the uniform electron 
gas is very far from magnetism in this density range. 
In solids near quantum critical points, the result is an 
overestimate of the magnetic moments and tendency to- 
ward magnetism (i.e. misplacement of the position of 
the critical point) due to neglect of the quantum critical 
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fluctuations. 36,37 

The present results for Na^CoC^ show a weak ferro- 
magnetic instability that is robust with respect to doping 
and structure (note the instability for the strained lat- 
tice). Based on this, and the experimentally observed 
renormalized paramagnetic state, it seems likely that 
Na^CoO^ is subject to strong ferromagnetic quantum 
fluctuations of this type, and that these are the reason 
for the disagreement between the LDA and experimental 
ground states. 

The effects of such quantum fluctuations can be de- 
scribed on a phenomenalogical level using a Ginzburg- 
Landau theory in which the magnetic properties defined 
by the LDA fixed spin moment curve are renormalized by 
averaging with an assumed (usually Gaussian) function 
describing the beyond LDA critical fluctuations. 38 Al- 
though a quantitative theory allowing extraction of this 
function from first principles calculations has yet to be 
established, one can make an estimate based on the LDA 
fixed spin moment curves as compared with experiment. 
In particular, Na^CoO^ shows a disagreement between 
the LDA moment and experiment equal to p = 1 — x, 
and has a very steeply rising LDA energy for moments 
larger than p. Thus one may estimate an r.m.s. am- 
plitude of the quantum fluctuations of £ w ap in fig, 
with 1/2 < a < 1, and most likely closer to 1. These 
are large values c.f. ZrZn 2 . It is therefore tempting to 
associate the superconductivity of Naa;Co02 • yH 2 with 
ferromagnetic quantum critical fluctuations. Considering 
the simple 2D Fermi surface, which consists of rounded 
hexagonal cylinders plus small sections, 28 and the ferro- 
magnetic fluctuations a triplet state like that originally 
discussed for Sr 2 RuC>4 (Rcf. 7) seems plausible. I now 
speculate about the ingredients in a spin fluctuation me- 
diated triplet superconducting state. 

Within a spin fluctuation induced pairing approach 
analogous to that employed for Sr 2 Ru04 the key ingre- 
dient is the integral over the Fermi surface of the k- 
dependent susceptibility with a function of the assumed 
triplet symmetry, 8,9 ' 17 i.e. in the simplest case, k-k'/fcfc'. 
For a Fermi surface in the shape of a circular cylin- 
der, radius kp, the needed integral is proportional to 
/j' d9cos(9)V(2kFsin(9/2)), where V(k) is the assumed 
pairing interaction. In any case, for a smooth variation 
of the spin fluctuations with k and a maximum at k~Q 
(ferromagnetic), the integral is roughly proportional to 
kp times the variation of V from k = to k = 2fcp. This 
latter variation depends on the detailed shape of V(k), 
but may be expected to cross over from being propor- 
tional to kp for small kp to proportional to kp for larger 
kp. Neglecting small Fermi surface sections, kp varies as 
p 1 / 2 . One possibility for V(k) is a function smoothly go- 
ing from a finite value at k = to near zero at the zone 
boundary (reflecting the rather weak antiferromagnetic 
instability relative to the ferromagnetic), with a size at 
k = given by the LDA ferromagnetic energy (oc p 2 ) 
or alternately a Hund's coupling (p independent) times £ 



(oc p). 

Within this p-wave scenario it would be quite interest- 
ing to measure the variation of the superconducting prop- 
erties of Na^CoO^ • J/H2O as a function of doping level. 
The above arguments imply a substantial model depen- 
dent variation up to the level where proximity to the crit- 
ical point suppresses T c , with the implication that still 
higher values of T c may be obtained. However, it should 
be emphasized that the mechanism and superconducting 
symmetry of Na^OoO^ -yH^O have yet to be established, 
and in fact, even conventional electron-phonon supercon- 
ductivity competing with spin fluctuations has not been 
excluded. 

I am grateful for helpful discussions with R. Asahi, 
I.I. Mazin, A.J. Millis, W.E. Pickett, S.S. Saxena and I. 
Terasaki. Some computations were performed using fa- 
cilities of the DoD HPCMO ASC and ARL centers. The 
DoD-AE code was used for some of this work. 
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TABLE I. LSDA spin magnetizations and energies for 
Naa;Co02. All quantities are on a per Co basis. Energies 
are in meV, spin moments are in fiB, FM denotes ferromag- 
netic and AF denotes the partially frustrated nearest neighbor 
AF configuration discussed in the text. M is the total spin 
magnetization, m is the magnetization inside the Co LAPW 
sphere, radius 1.95 Bohr. 
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FIG. 1. LDA fixed spin moment energy as a function of 
constrained spin magnetization of Na;cCo02 on a per Co atom 
basis for £=0.7 (x), £=0.5 (*) and a;=0.3 (+). The curves are 
spline interpolations as a guide to the eye. Note the breaks 
at 0.3, 0.5 and 0.7 \ib for £=0.7, 0.5 and 0.3, respectively. 
These correspond to the band gap between crystal field split 
Co d manifolds. The calculations were done keeping the struc- 
ture fixed and varying the Na site occupation via the virtual 
crystal method (see text). 
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